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Abstract 

All the well known exactly solvable s— wave potentials V{r), r G (0, oo) may be 
PT— symmetrically regularized (say, by a "small" imaginary shift of the coordinate 
axis) and re-defined as acting on the whole real line, r — > x—i e, x G {—oo, oo). There 
is no surprise that the new models are still exactly solvable and that their energies 
remain real. What is an unexpected result? The observed drastic and unpredictable 
changes in the forms of the new spectra themselves. 
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1 Introduction 



This lecture is mainly intended as a preliminary review and summary of the new ex- 
actly solvable models obtained by the author, step-by-step, during the recent months, 
and accessible already in the form of several short notes. 

In essence, the overall structure or "contents" of this text on non-Hermitian 
models might just parallel the page 297 of the review paper [1]. There, the exactly 
solvable Hermitian models have been listed and characterized by their so called shape 
invariance. 

The quoted list of potentials further decays in the two separate families. Let us 
call them C and for our present purposes. They are characterized by the respective 
Laguerre and Jacobi polynomial forms of their bound states. 

In the former family the first subcategory CX is defined on the whole axis, x G 
(—00,00). It contains just the quadratic oscillator V^'^'^\x) and the exponential 
Morse interaction V^^^\x). 

The second subcategory CR of the former family consists of the two best known 
three-dimensional solvable models, viz., of the Coulomb or Kepler V^^^\r) and the 
spherical harmonic V'^^^\r), both with r G (0, 00) of course. 

The latter family J' has been conveniently split in the three subcategories with 
two elements each. These subcategories are characterized by their ranges of coordi- 
nates. 

A particle "lives" on a finite interval in the first subcategory JY. Temporarily, 
let us skip and ignore this case completely. 

The second subcategory J'X is an opposite extreme with coordinates covering 
the whole real axis, x e IR. We may refer to its elements as Rosen Morse (x) 
and the scarf potential V'^^^^x). 

This agrees with the convention accepted in ref. [1]. Similarly, the respective 
Poschl- Teller and Eckart forces V'^^'^\x) and V'^^^\x) are simply determined as 
contained in the last, third subcategory J^R of the whole list. 

In what follows, I shall try to parallel this classification of the Hermitian solvable 
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models within the shghtly more general, non-Hermitian framework of the so called 
VT symmetric quantum mechanics [2] . In detail I shall only outline some of the most 
interesting observations, new results and open questions concerning these brand-new 
and still fairly unusual models. 

A unifying principle of all the present generalizations will be the condition of 
the so called VT symmetry [2]. This, in essence, means the commutativity of our 
Hamiltonians H with the product of the parity V and the operator T of time reversal. 
One has to note that in all the present time-independent models the latter operator 
symbolizes just a reflection T i — —iT oi the complex plane with respect to its real 
axis. 

2 Various types of the PT— symmetric potentials 

One of the first studies of the VT symmetric potentials [3] paid attention to the 
quartic (i.e., unsolvable) power-law forces and has been motivated by their close 
formal (and, first of all, Fourier-transformation) connection to their Hermitian power- 
law analogues. 

Similar questions can be, of course, much more easily studied within the domain of 
solvable interactions. Unfortunately, their explicit descriptions and/or classifications 
has only been initiated, quite recently, after a few parallel pioneering proposals of the 
VT symmetric quantum mechanics by Andrianov et al [4], Bender et al [2], Cannata 
et al [5] et al. 

2.1 Generalizations within the family C 
2.1.1 CX{PT) 

I shall omit here the detailed analysis of the VT symmetrized QO and EM models. 
My main reason is an idea of keeping this text as short as reasonable, supported 
by another important fact that the majority of the most relevant comments on the 
VT symmetrized category CX{VT) (and on its two models denoted, naturally, as 
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Y(ptqo)^^^ and V^^'^^^\x)) are already easily accessible in the respective publica- 
tions [6] and [7]. 

2.1.2 CR{VT) 

I have mixed reasons why 1 shall also skip the other two Laguerre-related models 
y{PTCK)^^-^ and V^™\x). 

The reason remains similar in the latter case only. Indeed, most of the facts 
I know about it did already appear in my letter [8]. In contrast to that, I (and, 
apparently, most of us) do not know virtually anything about the former model 
y^PTCK)!^^-^ yet. Only some indirect hints have appeared in a few loosely related 
publications [9]. I strongly believe that the situation will significantly improve soon 
[10]. 

2.2 Generalizations within the subcategory J X 

In the frame of our present denotation, the respective VT symmetrized models 
y{PTSP)(^^'^ and V''^'^^^\x) have been proposed in the recent short communica- 
tions [11] and [12]. Considerations based on the economy of space force me again to 
restrain from any re-copied formulae and repeated comments. 

Let me confirm this rule of conduct by one slightly discomforting exception. One 
can really feel puzzled by the latter innocent-looking, smooth and asymptotically 
vanishing potential which can be always uniformly bounded, \V^^'^^^\x)\ < 5. After 
its deeper and constructive study we were still not that much surprised by its ability 
of supporting a ground state in the weak coupling regime where, by our assumption, 
both its couplings were small of the order 0{5) at least. 

Nevertheless, one cannot help feeling caught our of guard by the observation that 
the related ground state energy may become also arbitrarily large, Eq — Oil/S"^) 
[12]. Similar paradoxes offer one of the reasons why the solvable models of the VT 
symmetric quantum mechanics [2] deserve a thorough and exhaustive analysis. 
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2.3 Generalizations within the subcategory J R 

Moving to a central point of this paper, let us point out that there exists one key 
difference between the previous oscillators and the "last two" PT symmetrized mod- 
els ^(^^^^)(r) and V^(^^^'^)(r) since both the Hermitian predecessors of these two 
forces are, generically, strongly singular in the origin. 

In the former model ^^^^^^^(r), its smoothness or regularization has been achieved 
by a small, local deformation of the integration path r = r(i), t & M near the origin. 
In addition, a complex rotation of one of the couplings by 7r/2 has been employed, 
in our preprint [13], as a guarantee of the VT symmetry. 

Such a complexification changes the spectrum and becomes a source of its several 
unusual features. A detailed inspection reveals, e.g., that an increase of the repulsion 
can lower (!) the energy. More details will be recollected here in sect. 3 below. 

In a less drastic approach to the singular forces we preserved the real values 
of the couplings in V^^^^^\r). We employed the "minimal" VT symmetrization 
r{t) — t — ie oi refs. [3] or [8]. Detailed properties of the resulting "last" solvable 
model have been described in the preprint [14]. At length they will also be discussed 
in sect. 4 below. 

Quotation marks in the word "last" mean that wc introduced one more model 
very recently [15]. The resulting force V^^^^^\x) of a Hulthen-type shape fits nicely 
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in the following scheme, 



symmetric 
straight — line 

YiPTSH)(^^^ 

in jCR [8] 



X r — —i exp i x 



symmetric 
straight — line 

in JR [14] 



I sinh r — —i exp i x 



periodic 

yiPTEM)^^^ 

on an arch 
in CX [7] 



periodic 

on an arch x{t), 
in JX [15] 



The vertical correspondence originates from the changes of variables and the hori- 
zontal arrows indicate the transition between the families C and J . One notices the 
similarities in the (symmetric or periodic) form of the functions V as well as in the 
(straight- line or bent-curve) shape of their domains r — r{t) or x — x{t) e C. 
More details will be provided in our final sect. 5. 



3 A generalization of the Eckart oscillator 

A deeper understanding of one-dimensional systems may be mediated by an analytic 
continuation of their real Schrodinger equation 

iP{r) = Eip{r), ^/'(±oo) = 0. (1) 

For one of the simplest particular models V{x) = a; A the loss of hermiticity at 
complex couplings proved more than compensated by the new insight in its solutions. 
E.g., its spectrum is given by a single multi-sheeted analytic function of A e C 
[16]. The same idea has been re- applied to the set of resonances in the cubic well 
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V{x) = LU x'^+Xx^ [17]. In the cubic case it was surprising to notice that the spectrum 
En{X) remains real for the purely imaginary couplings X — ig. The rigorous proof 
of this curious observation dates back to the late seventies [18]. It went virtually 
unnoticed for more than ten subsequent years. The phenomenon only re-entered 
the physical scene with Zinn- Justin and Bessis who, tentatively, attributed the strict 
absence of decay lmEn{i g) = to the mind-boggling real-symmetry-plus-imaginary- 
antisymmetry of the cubic force in question [19]. They also performed a number of 
numerical experiments, keeping in mind a paramount importance of this peculiar 
symmetry in field theory. There, it precisely coincides with the fundamental VT 
(i.e., parity-plus-time- reversal) invariance. According to Bender et al [20] this new 
type of symmetry might even replace the traditional requirement of hermiticity in 
many phenomenological models. 

Within the quantum mechanics itself the parallels between g and igx^ in- 
spired the numerical and semi-classical study of the generalized anharmonic forces 
V^'^\x) = LO + gx'^iixY with a variable real exponent 5 [6]. Within the related 
PT— symmetric branch of the "classical" quantum mechanics there appeared new 
perturbation series [21] and quasi-classical approximations [22], a new implementa- 
tion of supersymmetry [4] and the new types of spectra [23]. 

Among all the different models with the T'T— symmetrically broken parity one 
may distinguish, roughly speaking, between its "stronger" and "weaker" violation. 
The former group is formally characterized by the globally, asymptotically deformed 
paths of integration in eq. (1). An illustration may be provided by the elementary 
ground-state wave function = exp {—ix^-\- hx"^) of Bender and Boettcher [24] 
which ceases to be integrable on the real line of x. The integrability is only recovered 
after we bend both the semi-axes downwards, 

{x>l} {x = Qe-''P}, {x<-l} {x^-Qc''^'} 

with ^ ^ 1 and < </? < 7r/3. The wave function obviously corresponds to the 
quasi-exactly solvable potential V{x) = —9x^ — 12bix^ + ib'^x^ — Qix [24] and mimics 
the choice of 5 = 2 in the family V^^\x). The further growth oi d > 2 would make 
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both the asymptotical (^—wedges shrink and rotate more and more downwards in 
the complex plane. 

The second group of the VT symmetric examples with a "weaker" parity break- 
down does not leave the real axis of x at all (i.e., = 0, cf., e.g., [25]). This admits 
the more natural physical interpretation of the real physical coordinates. Such a form 
of the P— violation has been also implemented in several numerical and perturbative 
models. Their subclass which possesses elementary solutions is particularly instruc- 
tive since it incorporates all the so called shape invariant one-dimensional models of 
the ordinary quantum mechanics [12]. 

In both the groups of examples an overall VT symmetry of the Hamiltonian is, 
presumably, responsible for its real and discrete spectrum [6] . Cannata et al [5] were 
the first to notice that one of the various limits (5 — > oo of the power-law models 
with (/? — > 7r/2 — OiljS) becomes exactly solvable in terms of Bessel functions. This 
re-attracted attention to the related strongly deformed contours [26] . More recently, 
the same merit of an indirect formal parallel to the Hermitian square well has been 
also found for the standard real contour. The related 5 — > oo wave functions even 
degenerated to the Laguerre polynomials [7]. 

The latter unexpected emergence of the new exactly solvable model within the 
generalized, T'T— symmetric quantum mechanics encouraged our present study. In- 
deed, exactly solvable models are obviously best suited for analyses of methodical 
questions. In particular, the class of the PT— symmetrized shape invariant oscil- 
lators [12] does not seem to differ too much from its Hermitian counterpart. For 
an explicit analysis of the details of this correspondence one may simply recall the 
numerous explicit formulae available, e.g., in Table 4.1 of the review [1] or in the 
original factorization constructions [27] and their Lie-algebraic [28], operator [29] or 
supersymmetric [30] re-interpretations. 

Seemingly, one cannot expect any interesting new developments in the exactly 
solvable context. Fortunately, in the light of our recent remark on the spherical 
harmonic oscillator [8] non-trivial innovations may be expected in the domain of 
singular forces. Indeed, within the PT— symmetric quantum mechanics it is possible 
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to avoid some isolated singularities by a local deformation of the integration path. 
In particular, a strong repulsion in the origin (so popular in some phenomenological 
models [31] but fully impenetrable in one dimension) may be readily controlled by a 
suitable choice of the cut. 

In the present text we intend to re-attract attention to the singular forces. In 
eq. (1) we shall use the asymptotically real path of integration which is only locally 
deformed and non-Hermitian. We shall show that this innovative approach enables 
us to regularize the one-dimensional models via their suitable VT symmetrization. 

3.1 Eckart model 

Our particular attention will be paid to the exceptional s— wave potential 

sinh X sinhx 

with the strongly singular core. Usually attributed to Eckart [32], this model is 
solvable on the half-hne with x e (0, oo) and, conventionally, A> 1/2 and B > [1] . 
Its fixed value of the angular momentum £ = is in effect a non-locality which lowers 
its practical relevance in three dimensions. Here, we shall study its PT— symmetrized 
version with the purely imaginary coupling B = i(3. Besides the obvious relevance 
of such an exceptional complexified model with a strong singularity in quantum 
mechanics, an independent encouragement of our study is also provided by its obvious 
phenomenological and methodical appeal in the context of field theory, especially in 
connection with the so called Klauder phenomenon [33]. 

The local deformation of the integration path will enable us to forget about the 
strong singularity in the origin. This deformation will also admit the presence of the 
so called irregular components in 'il){x) ~ x^~^ near x — Q. They would be, of course, 
unphysical in the usual formalism [34]. 
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3.2 Solution revisited 

For all these reasons we have to re-analyze the whole Schrodinger equation anew. 
Our initial choice of the appropriate variables 

is still dictated by the arguments of Levai [30]. Then we insert V'^^'^''\x) in eq. (1) 
and our change of variables leads to its new form 

z{l - z) ip"{z) + [c - (a + 6 + 1)^] ip'{z) - ah cp{z) = (2) 

where 

c=l + 2w, a + b^2u + 2v + l, ab ^ (u + v){u + v + 1) + A(l - A) (3) 

and 

Av'^ ^2B-E, Au^ = -2B - E. (4) 

Our differential equation is of the Gauss hypergeometric type and its general solution 
is well known [35], 

^{z) = Ci • 2i^i(a, h- c; 2;) + C2 • z'^'^ 2^^1(0 + 1 - c, 6 + 1 - c; 2 - c; 2). (5) 

The first thing we notice is that our parameters a and h are merely functions of the 
sum u-\-v and vice versa, u-\-v — {a + h — l)/2. The immediate insertion then gives 
the rule (a — 6)^ = {2A — Vf and we may eliminate 

a^h±{2A-l). (6) 

We assume that our solutions obey the standard oscillation theorems [36] and become 
compatible with the boundary conditions in eq. (1) at a discrete set of energies, i.e., 
if and only if the infinite series 2-^1 terminate. Due to the complete a <-> 6 symmetry, 
we only have to distinguish between the two possible choices of C2 and Ci = 0. 

In the former case with the convenient h = —N (= non-positive integer) the 
resulting numbers a + b and u + v prove both real. Using the definition of B the 
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difference u — v = —ij3 /{u + v) comes out purely imaginary. The related terminating 
wave function series (5), i.e., 

is asymptotically normalizable if and only if m + v > 0. This condition fixes the sign 
in eq. (6) and gives the explicit values of all the necessary parameters, 

a^2A-N-l, u + v^A-N-l, u-v^-i- — ^- . (8) 

For all the non-negative integers N < N^ax < A — 1 the spectrum of energies is 
obtained in the following closed form, 

E^-\ («' + ^') =-(^-^-l)' + pr:^nj^, 7V = 0,1,...,7V^„.. (9) 

The normalizable wave functions become proportional to Jacobi polynomials, 

ip[z{x)] = canst. ■ P^^''^/'^(cothx). (10) 

Before we start a more thorough discussion of this result we have shortly to return 
to the second option with Ci = in eq. (5). Curiously enough, this does not 
bring us anything new. Although the second Gauss series terminates at the different 
b — c — 1 — N , the factor z^~'^ changes the asymptotics and one only reproduces 
the former solution. All the differences prove purely formal. In the language of our 
formulae one just replaces u by —u in (and only in) both equations (7) and (8). No 
change occurs in the polynomial (10). 

3.3 Spectrum 

The new spectrum of energies seems phenomenologically appealing. The separate 
N—ih. energy remains negative if and only if the imaginary coupling stays sufficiently 
weak, 13'^ < [A — N — 1)^. Vice versa, the highest energies may become positive, 
with E = E{Nmax) growing extremely quickly whenever the value of the coupling 
A approaches its integer lower estimate 1 + N^ax from above. In this way, even 
a weak VT symmetric force V^^^''\x) is able to produce a high- lying normahzable 
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excitation. This feature does not seem connected to the presence of the singularity 
as it closely parallels the similar phenomenon observed for the symmetric Rosen- 
Morse oscillator which remains regular in the origin [12]. Also, in a way resembling 
harmonic oscillators the distance of levels in our model is safely bounded from below. 
Abbreviating D — A — N — 1 — A^fff.ctive > its easy estimate 

(useful, say, in perturbative considerations) may readily be improved to Ej^ — E^_i > 
at small D ^ 1, to En — E^-i > '2D at large D 1 and, in general, to an 
algebraic precise estimate obtainable, say, via MAPLE [37]. 

Let us emphasize in the conclusion that the formulae we obtained are completely 
different from the usual Hermitian s— wave results as derived, say, by Levai [30]. He 
has to start from the regularity in the origin which implies an opposite sign in eq. (6). 
This must end up with the constraint B > 0. Moreover, the size of B would hmit 
the number of bound states. In the present VT symmetric setting, a few paradoxes 
emerge in this comparison. Some of them may be directly related to the repulsive 
real core in our V^^'^'^\x) with imaginary B. Thus, one may notice that the increase 
of the real repulsion lowers the A^— th energy. In connection with that, the number 
of levels grows with the increase of coupling A. In effect, the new bound-state levels 
emerge as decreasing from the positive infinity (!). At the same time, the presence 
of the imaginary B — i(3 shifts the whole spectrum upwards precisely in the manner 
known from the non-singular models. 

4 A generalization of the Poschl- Teller potential 

Among all the exactly solvable models in quantum mechanics the one-dimensional 
Schrodinger equation (1) with one of the most elementary bell-shaped potentials 
V^^^\r) — G/cosh^r is particularly useful. Its apphcations range from the analyses 
of stability and quantization of solitons [38] to phenomenological studies in atomic 
and molecular physics [39], chemistry [40], biophysics [41] and astrophysics [42]. Its 
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appeal involves the solvability by different methods [30] as well as a remarkable 
role in the scattering [43]. Its bound-state wave functions represented by Jacobi 
polynomials offer one of the most elementary illustrations of properties of the so 
called shape invariant systems [1] . The force V'^^^^ (r) is encountered in the so called 
VT symmetric quantum mechanics [2] where it appears as a Hermitian super-partner 
of a complex "scarf model [11]. 

Curiously enough, it is not too difficult to extend the exact solvability of the 
potential V^^^\r) to all its "spiked" (often called Poschl-Teller [44]) shape invariant 
generalizations 

Vi^Hr) = -^^ + ?^. (11, 
cosh r smh r 

In a way resembling the preceding section, our new one-dimensional Schrodinger eq. 
(1) is also too singular at B{B — 1) 7^ 0. The force ^^^^^(r) must be confined to the 
semi-axis, r e (0, cxo). This makes the "improved" Poschl-Teller model (11) much 
less useful in practice since its higher partial waves are not solvable. 

4.1 Regularization 

We may repeat that the impossibility of using eq. (11) in three dimensions (or on 
the whole axis in one dimension at least) is felt unfortunate because the singular 
potentials themselves are frequently needed in methodical considerations [33] and in 
perturbation theory [45] . They are encountered in phenomenological models [46] and 
in explicit computations [31] but not too many of them are solvable [47]. 

We feel mainly inspired by the pioneering letter [6] where Bender and Boettcher 
modified the harmonic oscillator V^^'^\r) = by a complex downward shift of its 
axis of coordinates, 

r — x — ie, xe(— 00,00). (12) 

The VT symmetry of their model V^(^^) (x) = \/(^'^) (x — ic) — — 2icx — (? means 
its invariance with respect to the simultaneous refiection x — > —x and complex 
conjugation i — > —i. Various other complex interactions have been subsequently 
tested and studied within this framework. 
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Our study [8] of the three-dimensional VT symmetric harmonic oscillator offers 
the details of our present key idea. The shift (12) has been employed there as 
a source of a regularization of the strongly singular centrifugal term. As long as 
l/{x — ief = (x + isY /{x^ + e^Y at any £ 7^ 0, this term remains nicely bounded 
in a way which is uniform with respect to x. Without any difficulties one may work 
with V^^^^\x) = r^{x)+£{£+l)/r^ (x) on the whole real line of x. In what follows 
the same idea will be applied to the regularized Poschl-Teller-like potential 

= - te), 0<e< 7i/2. 

This potential is a simple function of the Levai's [30] variable g{r) = cosh2r. As 
long as g{x — i e) = cosh 2x cos 2e — i sinh 2x sin 2e, the new force is VT symmetric 
on the real line of x e (—00, 00), 

Due to the estimates | sinh^ (x — ie)]'^ — sinh^ x cos^ e + cosh^ xsin^ e — sinh^ x + sin^ e 
and I cosh^(a; — is)\'^ = sinh^x + cos^e the regularity of V'^^^'^\x) is guaranteed for 
all its parameters e G (0,7r/2). 

4.2 Solutions 

In a way paralleling the three-dimensional oscillator the mere analytic continuation 
of the s— wave bound states does not give the complete solution. One must return to 
the original differential equation (1). There we may conveniently fixA-|-l/2 = Q;>0 
and S - 1/2 = /3 > and write 

/ (3^-1/4 a^-l/4\ , ^ , 

\-^+ .2 . txTtW'^^)^^^^^)' r{x)^x-ie. (13) 

y dx'^ smh r{x) cosh r{x) ) 

This is the Gauss differential equation 

z{\ + z) ip"{z) + [c+{a + b+ l)z] ip'{z) + ah ^{z) = (14) 
in the new variables 

'il){x) ^ z^{l + zy(f{z), 2; = sinh^r(a;) 
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using the suitable re-parameterizations 



a' = {2v - l/2)^ 



= (2/x - l/2)^ 



2//+l/2 = c, 2// + 2i/ = a + &, E^-{a-hf. 



In the new notation we have the wave functions 



^(x) = sinh"^+^/2[r(a;)] cosh'^+^Z^^^^^)] ^[z{x)] 



(15) 



with the sign ambiguities r = ±1 and u = ±1 in 2/x = r/? + 1/2 and 2v — aa-\- 1/2. 
This formula contains the general solution of hypergeometric eq. (14), 

^{z) = Ci 2i^i(a, b- c; -z) + C2Z^'' 2Fi{a + 1 - c,h + I - c-2 - c- -z). (16) 

The solution should obey the complex version of the Sturm-Liouville oscillation the- 
orem [36]. In the case of the discrete spectra this means that we have to demand 
the termination of our infinite hypergeometric series. This suppresses an asymptotic 
growth of 'il){x) at X ^ ±oo. 

In a deeper analysis let us first put C2 = 0. We may satisfy the termination 
condition by the non-positive integer choice of 6 = —N. This implies that a — 
N + 1 + aa + Tf3 is real and that our wave function may be made asymptotically 
(exponentially) vanishing under certain conditions. Inspection of the formula (15) 
recovers that the boundary condition '0(±oo) = will be satisfied if and only if 



The closed Jacobi polynomial representation of the wave functions follows easily, 



Now we have to return to eq. (16) once more. A careful analysis of the other 
possibility Ci = does not recover anything new. The same solution is obtained. 



1 < 2Ar + 1 < 2Nmax + 1 < -(^ci - t(3. 




The final insertions of parameters define the spectrum of energies. 



E = -(2A^ +l + aa + T(5f < 0. 



(17) 
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with r replaced by — r. Wc may keep C2 = and mark the two independent solutions 
by the sign r. Once we define the maximal integers Njj^^J which are compatible with 
the inequality 

2Ar(,^;) + l<-aa-T/3 (18) 
we get the constraint N < . The set of our main quantum numbers is finite. 

4.3 Paradoxes 

Let us now compare our final result (17) with the known £ = formulae for s 
waves [30] . An additional physical boundary condition must be imposed in the latter 
singular limit [48]. This condition fixes the unique pair a — —1 and r = +1. Thus, 
the set of the s— wave energy levels Ej^ is not empty if and only if a — /3 > 1. In 
contrast, all our e > potentials acquire a uniform bound \ V^^^'^{x)\ < const < 00. 
Due to their regularity, no additional constraint is needed. Our new spectrum E^''^^ 
becomes richer. For the sufficiently strong couplings it proves composed of the three 
separate parts, 

Ep-^ < 0, 0<N< TVK-), a + (3>l, 

£^{7'^U0, 0<N<N^+\ a>P + l, (19) 

Ei^'~^ < 0, < iV < N^^-\ P>a + 1. 

The former one is non-empty at A + B > 1 (with our above separate conventions 
A > —1/2 and B > 1/2). Concerning the latter two alternative sets, they may 
exist either at A > B or at B > A + 2, respectively. We may summarize that in a 
parallel to the VT symmetrized harmonic oscillator of ref. [8] we have the N^^^^ + 1 
quasi-odd or "perturbed" , analytically continued s— wave states (with a nodal zero 
near the origin) complemented by certain additional solutions. 

In the first failure of a complete analogy the number N^^~^ -|- 1 of our quasi- 
even states proves systematically higher than N^^^^^ + 1, especially at the larger 
"repulsion" 1. This is a certain paradox, strengthened by the existence of 

another quasi-odd family which behaves very non-perturbatively. Its members (with 
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the ground state 'ijjQ^'^\x) = cosh'^^^[r(x)] sinh^~^[r(a;)] etc) do not seem to have any 
s— wave analogue. They are formed at the prevalent repulsion B > A+2 which is even 
more counter- intuitive. The exact solvability of our example enables us to understand 
this apparent paradox clearly. In a way characteristic for many VT symmetric 
systems some of the states are bound by an antisymmetric imaginary well. The whole 
history of the VT symmetric models starts from the purely imaginary cubic force 
[19] after all. A successful description of its perturbative forms V{x) — lox^ + i\x^ 
is not so enigmatic [18] , especially due to its analogies with the real and symmetric 
V{x) — ojx^^Xx'^ [17]. The similar mechanism creates the states with (cr, r) = (+, — ) 
in the present example. A significant novelty of our new model V^^^"^^ {x) lies in the 
dominance of its imaginary component at the short distances, a; ~ 0. Indeed, we may 
expand our force to the first order in the small £ > 0. This gives the approximation 

1 smh^ix + ie) 1 ^. coshx n, 

27 ^ = 7 = 2— + 2«£ 3- + 0(£^). (20) 

sinh {x — is) (sinh x + sin eY sinh x sinh x 

We see immediately the clear prevalence of the imaginary part at the short distances, 
especially at all the negligible A — 0{e^). 

An alternative approach to the above paradox may be mediated by a sudden 
transition from the domain of a small £ ~ to the opposite extreme with e 7r/2. 
This is a shift which changes cosh a; into sinhx and vice versa. It intertwines the 
role of a and /? as a strength of the smooth attraction and of the singular repulsion, 
respectively. The perturbative/non-perturbative interpretation of both our quasi- 
odd subsets of states becomes mutually interchanged near both the extremes of e. 

The dominant part (20) of our present model leaves its asymptotics comparatively 
irrelevant. In contrast to many other PT symmetric models as available in the 
current literature our potential vanishes asymptotically. 

An introduction and analysis of continuous spectra in the VT symmetric quantum 
mechanics seems rendered possible at positive energies. This question will be left 
open here. In the same spirit of a concluding remark we may also touch the problem 
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of the possible breakdown of the VT symmetry. This has recently been studied on the 
background of the supersymmetric quantum mechanics [4]. In our present solvable 
example the violation of the VT symmetry is easily mimicked by the complex choice 
of the couplings a and /3. Due to our closed formulae the energies will still stay real, 
provided only that Im [aa + r/3) = 0. 

5 VT symmetric Hulthen potential 

Quantum mechanics is often forced to formulate its predictions numerically. Its ex- 
actly solvable models are scarce. Only their subclass in one dimension is broader 
and, in this sense, privileged and exceptional. It involves the harmonic oscillator and 
Morse potentials (with bound states expressible in terms of the Laguerre polynomi- 
als) as well as several other, less popular models solvable in terms of the polynomials 
of Jacobi etc (cf., e.g., review [1] for more details). 

One encounters a high degree of analogy between the real and complex forces for 
the exactly solvable models based on the use of Jacobi polynomials. An appropri- 
ate exactly solvable VT symmetrization of both the one-dimensional shape- invariant 
models with the "canonical" Rosen-Morse and scarf-hyperbohc shapes has been re- 
cently proposed and analyzed in refs. [11]. 

As we already noticed, a word of warning comes from the phenomenologically 
most appealing Coulomb interaction. Up to now, its only available VT symmetric 
simulation proves merely partially solvable [9]. In this section we intend to offer 
a partial remedy. We shall derive and describe a complete and exact solution for 
an appropriate VT symmetric complexification of the singular phenomenological 
Hulthen potential [43] which is known to mimic very well the shape of the Coulombic 
force in the vicinity of its singularity. 
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5.1 Method: The Liouvillean change of variables 



In the first step let us recollect that in the spirit of the old Liouville's paper [49] the 
change of the (real) coordinates (say, r ^) in Schrodinger equation 

^2 



X{r) = x(r) 



(21) 



may sometimes mediate a transition between two different potentials. It is easy to 
show [50] that once we forget about boundary conditions one simply has to demand 
the existence of the invertible function r = r(^) and its few derivatives r'(^), r"(^), . . . 

in order to get the explicit correspondence between the two bound state problems, 
viz., original (21) and the new Schrodinger equation with the known wave functions 



(22) 



generated by the new interaction 



^(0 -E^ [r'iO]" {WHO] + + l 





' 1 


\r"'{0] 


[r'iOl 


2 





(23) 



One might recall the well known mapping between the Morse and harmonic Laguerre- 
related oscillators as one of the best known exphcit illustrations. For it, the necessary 
preservation of the correct physical boundary conditions is very straightforward to 
check [43]. 

In the Jacobi-polynomial context the Liouvillean changes of variables have been 
applied systematically to all the Hermitian models (cf. Figure 5.1 in the review [1] 
or ref. [51] for a more detailed illustration). A similar thorough study is still missing 
for the VT symmetric models within the same class. 

In the present letter we shall try to fill the gap. For the sake of brevity we shall 
only restrict our attention to the VT symmetric initial eq. (21) with the Poschl- Teller 
potential studied and solved exactly in our recent preprint [14], 



W{r) = 



(3^ - 1/4 



a 



1/4 



r — X — le, 



sinh^ r cosh^ r ' 
The normalizable solutions are proportional to the Jacobi polynomials. 



CXD, OO 



(24) 



X(r) ^sinh^'^+^/Vcosh^^+^/V P^^'^''"") (cosh 2r) 
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at all the negative energies —k'^<0 such that 

K = kI^'^^ = -aa - T/3 - 2n - 1 > 0. 

These bound states are numbered by n = 0,1, . . . ,n^'^J and by the generalized 
parities a = ±1 and r = ±1. 

We may note that our initial VT symmetric model (21) remains manifestly regu- 
lar provided only that its constant downward shift of the coordinates r = r(j;) — x—ie 
remains constrained to a finite interval, e e (0, 7r/2). In this sense our initial model 
(24) is closely similar to the shifted harmonic oscillator. At the same time, one still 
misses an analogue of a transition to its Morse-like final partner V{^). In a key step 
of its present construction let us first pick up the following specific change of the axis 
of coordinates, 

sinhr(j;)(^) = —ie^^, ^ — v — iu. (25) 

The main motivation of such a tentative assignment lies in the related shift and 
removal of the singularity (sitting at r = 0) to infinity {u —>■ +oo). In fact, one 
cannot proceed sufficiently easily in an opposite direction, i.e., from a choice of a 
realistic V{^) to a re-constructed r(^). This is due to the definition (23) containing 
the third derivatives and, hence, too complicated to solve. 

We shall see below that we are quite lucky with our purely trial and error choice 
of eq. (25). Firstly, we already clearly see that the real line of x becomes mapped 
upon a manifestly VT symmetric curve ^ — v — iu in accordance with the compact 
and invertible trigonometric rules 

sinh X cos e — sinv, 
cosh X sin £ = e" cos v, 

i.e., in such a way that 

V = arctan {^) = v^,) e [v^-^),v^^{) = (-f + £, f - , 
u — ui^x) = i 111 (^sinh^ X -\- sin^ . 

This relationship is not too different from the Morse-harmonic equivalence studied 
in ref. [7] . Our present path of ^ is a very similar down-bent arch which starts in its 
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left imaginary minus infinity, ends in its right imaginary minus infinity while its top 
lies at X — V — and —u — — ti(o) = In 1/ sine > 0. The top may move towards the 
singularity in a way mimicked by the diminishing shift £ — > 0. Indeed, although the 
singularity originally occurred at the finite value r ^ 0, it has now been removed 
upwards, i.e., in the direction of —u — > +oo. 

5.2 Consequences 

The first consequence of our particular change of variables (25) is that it does not 
change the asymptotics of the wave functions. As long as r'(^) = itanhr(^) the 
transition from eq. (21) to (1) introduces just an inessential phase factor in ^(^). 
This implies that the normalizability (at a physical energy) as well as its violations 
(off the discrete spectrum) are both in a one-to-one correspondence. 

The explicit relation between the old and new energies and couplings is not too 
complicated. Patient computations reveal its closed form. With a bit of luck, the 
solution proves non-numerical. The new form of the potential and of its binding 
energies is derived by the mere insertion in eq. (23), 

^k) = (t4«f + t^. ^ = ''^- PS) 

At the imaginary ^ and vanishing ^4 = this interaction coincides with the Hulthen 
potential. 

In the new formula for the energies one has to notice their positivity. This is 
extremely interesting since the potential itself is asymptotically vanishing at both 
ends of its integration path. One may immediately recollect that a similar paradox 
has already been observed in a few other VT symmetric models with an asymptotic 
decrease of the potential to minus infinity [24, 52] . 

The exact solvability of our modified Hulthen potential is not yet guaranteed at 
all. A critical point is that the new couplings depend on the old energies and, hence, 
on the discrete quantum numbers n, a and r in principle. This could induce an 
undesirable state-dependence into our new potential. Vice versa, the closed solvabil- 
ity of the constraint which forbids this state-dependence will be equivalent to the 
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solvability at last. 

A removal of the obstacle means in effect a transfer of the state-dependence (i.e., 
of the n—, a— and r— dependence) in 

A = A{a) = 1 - C (= A + B) = _ ^2 

from C to j3. To this end, employing the known explicit form of k we may re- write 

C = C{a, T, n) = {aa + 2n + l){aa + 2n + 1 + 2t(5). (27) 

This formula is linear in Tj3 and, hence, its inversion is easy and defines the desirable 
state-dependent quantity j3 = (3{a,T,n) as an elementary function of the constant 
C. The new energy spectrum acquires the closed form 



E = E{a, T,n) ^ A + B + - 



1 r „ . A + B 



aa + 2n + 1 



(28) 



aa + 2n + l. 

Our construction is complete. The range of the quantum numbers n, a and r remain 
the same as above. 

In the light of our new result we may now split the whole family of the exactly 
solvable VT symmetric models in the two distinct categories. The first one "fives" on 
the real line and may be represented or illustrated not only by the popular Laguerre- 
solvable harmonic oscillator [8] but also by our initial Poschl- Teller Jacobi-solvable 
force (24). 

The second category requires a narrow arch-shaped path of integration which all 
lies confined within a vertical strip. It contains again both the Laguerre and Jacobi 
solutions. The former ones may be represented by the complex Morse model of ref. 
[7]. Our present new Hulthen example offers its Jacobi solvable counterpart. The 
scheme becomes, in a way, complete. 

The less formal difference between the two categories may be also sought in their 
immediate physical relevance. Applications of the former class may be facilitated by 
a limiting transition which is able to return them back on the usual real line. In 
contrast, the second category may rather find its most useful place in the methodical 
considerations concerning, e.g., field theories and parity breaking [53]. Within the 
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quantum mechanics itself an alternative approach to the second category might also 
parallel studies [5] of the "smoothed" square wells in non-Hermitian setting. 

In the conclusion let us recollect that the VT symmetry of a Hamiltonian replaces 
and, in a way, generalizes its usual hermiticity. This is the main reason why there 
existed a space for a new solvable model among the singular interactions. An exactly 
solvable example with an "intermediate" (i.e., hyperbola-shaped) arc of coordinates 
remains still to be discovered. Indeed, this type of a deformed contour has only been 
encountered in the "quasi-solvable" (i.e., partially numerical) model of ref. [24]) 
and in the general unsolvable forces studied by several authors by means of the 
perturbative [18], numerical [54] and WKB [55] approximative techniques. 
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